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The Perfect Laplace Operator for Non-Trivial Boundaries 

S. Hauswirth a 

institute of Theoretical Physics, University of Bern 
Sidlerstrasse 5, CH-3012 Bern, Switzerland 

The application of Renormalization Group (RG) methods to find perfect discretizations of partial differential 
equations is a promising but little investigated approach. We calculate the classically perfect fixed-point Laplace 
operator for boundaries of non-trivial shape analytically and numerically and present a parametrization that can 
be used for solving the Poisson equation. 



1. INTRODUCTION 

In Lattice QCD, the perfect action provides a 
discrete formulation of the theory which is free 
from lattice artifacts jjj. Few attempts have so 
far been made to examine perfect discretizations 
also for classical problems like the numerical so- 
lution of differential equations J|-||]. We present 
an exploratory study of a perfect differential op- 
erator in presence of non-trivial boundary con- 
diations, as they appear in classical field theory 
applications. This work is in some more detail 
also presented in S. 

Our toy model is the Laplace equation on 
a two-dimensional region with fixed boundaries 
which are built from vertical and horizontal lines 
as shown in Fig. |l|. 



- the classically perfect fixed point operator - is 
completely free of discretization errors. Perfect 
operators have non-zero couplings spreading over 
an infinite volume, so for practical use they have 
to be truncated. Suppose we want to use an oper- 
ator with only nearest- and next-to-nearest neigh- 
bour couplings. Then we have to provide a set of 
operators which differ in the position of the cen- 
ter point with respect to the boundary (e.g. for 
point A in Fig. ^ we need another operator than 
for point B) as symbolized by the icons in Fig. |^: 
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Figure 2. The set of Laplace operators Pi(ri,r 2 ) 
near straight boundaries, including couplings 
with r t < 2 (i = 1,2). 



Figure 1. A two-dimensional bounded region. 

Solving such a problem numerically with an 
improved discretization of the Laplace operator 
reduces the need for high resolution of the un- 
derlying mesh. The most improved discretization 



We construct such a set of operators with the 
perfect action formalism. In order to check how 
well this construction works, we examine how 
the boundary influences the perfect Laplacian. 
We find that not only the operator itself is lo- 
cal (i.e. its couplings decrease exponentially with 
distance), but also the effect of the boundary on 
the operator (so that for example in point B we 
don't have to care about any boundaries). 
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2. THE PERFECT LAPLACIAN 



3. SQUARE-SHAPED BOUNDARIES 



The perfect Laplace operator in d dimensions 
can be calculated from the fixed point action of a 
free real scalar field with the continuum action 



A c 
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d d x d ll (j)(x)d l _ l (l)(x). 



(1) 



The corresponding equation of motion is the 
Laplace equation. A general discretization of the 
action (Q) contains terms which couple the field 
at one lattice site to the field at another one: 



A{4>) 



(2) 



Consider a square region with 4> — on the 
boundaries and side length Na. In this case, the 
perfect Laplacian can be calculated analytically. 
The easiest way to construct it in d = 2 is from 
symmetries: Extend the field cf> onto the whole 
plane by periodically mirroring it at the boundary 
with alternating sign 

<f>(x 1 ,x 2 ) = (-l) k i+ k ><f>((- 



-1) 



21N), 



(5) 



for any hi = 0, 1 and h G Z, (i = 1, 2). In Fig. g, 
the square in the center shows the region we are 
interested in, the field 4> is mirrored at all straight 
lines and its sign is given explicitly. 



with the couplings p(r), r = . . . , r^). For the 
standard Laplacian in two dimensions, p(0) = 4 
and p(r) = —1 (V|r| = 1). A RG transformation 
of the form 



n 



(3) 



relates the actions A((f>) on a fine and A' (x) on a 
coarse lattice with lattice units a and 2a, respec- 
tively. C is a normalization constant and T the 
blocking kernel 



(4) 



where k is a parameter used for optimizing the 
fixed point Laplacian in terms of locality. We 
use the kernel (||) with k — 2 in the following, 
as this gives a very local operator (the couplings 
decrease exponentially cx e~ jr with a large decay 
coefficient 7 « 3.5) [6,0|. The fixed point of the 
RG transformation (3) can be calculated analyt- 
ically l|. In momentum space, the result is 
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The couplings p* (n) of the fixed point Laplacian 
are found by Fourier transformation. From this 
result we derive the element p\ of the set of op- 
erators in Fig. H 



Figure 3. The field is periodically mirrored on the 
walls of squares. 



In the perfect Laplace equation for the whole 
plane (with the perfect Laplacian p*(n)), the sum 
over all lattice points splits up into sums over A^- 
squares 



'■>n+r 



1 AT-1 

EEE 

h£Z k i= Si=0 



p*{r-n)cp r =0,(6) 



where the variable r running over the lattice 
points is given by n = (— l) ki Si — ki + 2kN for i = 
1,2. Using the symmetry p*[—r%, r 2 ) =p*(ri,r2) 
and the relation (||), the perfect Laplacian for 
square boundaries is given by the infinite sum of 
couplings 

1 

p*(n,r) = £(-l) fel+fc V*(r+2fcn+fc+2Z7V), 
Uez k i= o 

with both n and n + r lying inside the bound- 
ary. The sum converges very fast because of the 
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Figure 4. The couplings p*(ni,n' 2 ) of the perfect 
Laplacian near a wall on a logarithmic scale with 
preserved sign. 



locality of p*{r). From this result we derive the 
operators P2~pe- 

The perfect Laplace operator for square bound- 
aries can also be recieved by performing the RG 
transformation analytically in Fourier space, us- 
ing basis functions V^q which enforce the condi- 
tion 4> — on the boundary. The result is p 
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The couplings of the perfect Laplacian are found 
by reverse Fourier transformation p*(n,n') = 
i/^ 2 S Q *g(»)*QK) ' P*(Q) with * Q (n) = 
CQi("l) • ZqiM and 



P (r, \ = I v / 2sin(Q i (n i + \)) ,Qi^n, 



(7) 



(Qi = kiw/N, ki = 1,...,N.) A compari- 
son of the results for the two different construc- 
tions shows that they provide the same couplings. 
They are also in agreement with . 



4. CONCAVE CORNERS 

In the case where the boundaries form a con- 
cave corner (as in P7-P13), the perfect operator 
can't be calculated analytically. Therefore we 
perform a numerical RG transformation itera- 
tively on a N = 32 lattice with L-shaped bound- 
aries as shown in Fig. |^: 
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Figure 5. The lattice used to calculate the perfect 
Laplacian near a concave corner. Points denote 
the fine and boxes the coarse lattice sites. 



For quadratic actions and blocking kernels, the 
RG transformation (||) can be written as a mini- 
mizing condition for the fine field <f>: 

A'( X ) = mm [Atf) + T ( X , </>)] + const. (8) 
<p 

To get close to a fixed point, we have to iterate 
this RGT step. The results of the previous step 
- which are the couplings of the resulting coarse 
action - are then used as an input for the next 
step, that is as a new starting guess for the fine 
action A(4>). After O(20) iterations, we find a 
very close approximation to the couplings of the 
fixed point action. 

We look for the couplings of the coarse field 
inside the large dashed box. Outside the small 
dashed box, we use in every RGT step the oper- 
ators pi~pe for the fine action. For the couplings 
inside the small dashed box, we use the standard 
Laplacian in the first iteration, and afterwards 
the result of the previous iteration. As a result 
we get the operators P7-/O13, completing our set. 
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5. RESULTS 

The effect of the boundaries on the fixed-point 
Laplacian is indeed highly local. This can be seen 
from the symmetry construction described above: 
In the even simpler case of a single wall along the 
second axis the actual operator p*(n,r) in terms 
of the translation invariant operator p*(r) (i.e. 
the operator in absence of boundaries) is given 
by 



p*(n,r) = p*(r) - p*(n + 2ni + l,r 2 ), 



(9) 



where ni + 1/2 denotes the distance from the wall 
in lattice units (ni=0,l,. . . ). From the argument 
of the last term we see that the effect of the wall 
decreases twice as fast with distance as the cou- 
plings of p*(r). The size of this difference term 
in dependence of the distance from the wall is 
shown in Fig. [(]. For n± > 3 it is below machine 
precision. 
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Figure 6. The relative difference Ap* = 
log | (p* (n, r) — p* (r)) / p* (r) | between the cou- 
plings near a wall and the respective couplings 
of the perfect Laplacian on an infinite lattice. 



When truncating the operators to nearest- and 
next-to-nearest neighbour couplings, one has to 
impose some conditions on the remaining cou- 
plings in order to ensure basic properties of the 
operator. These conditions are the correct disper- 
sion relation E(k) — \p\ and the limit for small 



momenta Ym\ q ^o p{q) — q 2 , leading to the sum 

rules En,r 2 / 5 ( ri > r 2) = and E n ,r 2 ( r l + r l) ■ 
p(ri,r 2 ) = -4. 

To check the quality of our parametriza- 
tion of the perfect Laplace operator, we nu- 
merically solved the lattice Poisson equation 
A<j) = f for square boundaries with the stan- 
dard (4, — 1, — 1, — 1)-Laplacian and with our 
parametrized Laplace operator made up from pi- 
pe. For the resulting field configuration, we de- 
termined the potential energy E = f J2 n 
for the field configuration </j( s * d ) computed with 
the standard Laplacian and E — J2 n 4>n V ^ + 
F 2 /3kV for the configuration </>( pp ) computed 
with the parametrized perfect Laplacian. The 
additional term F 2 /3kV in the potential energy 
comes from the constant C in the RG transfor- 
mation (||) for non-zero source § . 

Fig. shows a plot of the relative error AE = 
(S-.E;^ ™* ))/^^ ™* ) as a function of the inverse 
lattice volume 1/V = 1/N 2 . The error for our 
parametrization is proportional to 1 /N 2 and is for 
any lattice size smaller by a factor of about 180 
than the error of the standard Laplacian. This 
factor can probably still be increased by tuning 
the normalization procedure of the truncated cou- 
plings. 
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Figure 7. Relative error for the numerical solution 
of our test problem. 
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6. THE POISSON EQUATION 
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equation of motion 



FP 



where the fixed point source is given by 



J FP (k)=^(k)^-J(k). 



Q{k) ~, 



(12) 



(13) 



In the perfect Poisson equation (|12|) appears 
again the perfect Laplacian p*(r) which is de- 
scribed by our parametrization. 
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7. CONCLUSION 

The very basic examinations presented here 
show that the fixed-point approach is also ap- 
plicable to classical field theory problems. We 
show that the fixed point Laplace operator does 
not care about boundaries which are a few lat- 
tice spacings away, and therefore it is possible to 
give a simple parametrization which works well in 
applications. Whether this approach is also fruit- 
ful for more complicated problems remains to be 
investigated. 
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